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We study a pure gluon plasma in the context of quasiparticle models, where the plasma is con- 
sidered as an ideal gas of massive bosons. In order to reproduce SU(3) gauge field lattice data 
within such a framework, we review briefly the necessity to use a temperature-dependent gluon 
mass which accounts for color interactions between the gluons near Tc and agrees with perturba- 
tive QCD at large temperatures. Consequently, we discuss the thermodynamics of systems with 
temperature-dependent Hamiltonians and clarify the situation about the possible solutions pro- 
posed in the literature to treat consistently those systems. We then focus our attention to two 
possible formulations which are thermodynamically consistent and we extract the gluon mass from 
the equation of state obtained in SU(3) lattice QCD. We find that the thermal gluon mass is sim- 
ilar in both statistical formalisms. Finally, an interpretation of the gluon plasma as an ideal gas 
made of glueballs and gluons is also presented. The glueball mass is consistently computed within a 
relativistic formalism using a potential obtained from lattice QCD. We find that the gluon plasma 
might be a glueball-rich medium for T < 1.13 Tc and suggest that glueballs could be detected in 
future experiments dedicated to the quark-gluon plasma. 
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I. INTRODUCTION 

A. Generalities and lattice data 

It is expected that, at high enough temperatures or 
densities, a phase transition from hadronic matter to 
quark-gluon plasma will occur. As early as 1975, Collins 
and Perry suggested that the dense nuclear matter at the 
center of neutron stars could consist in deconfined quarks 
and gluons [ij. In 1980, Shuryak studied the nuclear mat- 
ter at high temperatures and introduced the terminology 
"quark-gluon plasma" in analogy with similar phenom- 
ena in atomic physics Qj- Beside its intrinsic interest, 
knowing the equation of state of a quark-gluon plasma is 
needed to predict the evolution of stars and to know for 
example if a neutron star can go through a quark phase or 
just collapse into a black hole. This knowledge is also im- 
portant to predict when our Universe hadronized, since it 
is also believed that it was a quark-gluon plasma within a 
few /is after the Big-Bang. The hadronic matter/quark- 
giuon plasma phase transition, predicted by Quantum 
Chromodynamics (QCD), is studied experimentally at 
the Relativistic Heavy Ion Collider (RHIC) ^ and will 
be also studied in the future at the Large Hadron Collider 
(LHC). The results obtained at RHIC so far suggest that 
the quark-gluon plasma behaves like an almost perfect 
fluid instead of a weakly interacting gas, indicating that 
interactions are still quite large after the phase transi- 
tion, for T > Tc {Tc denotes the critical temperature 
of QCD). From a field-theoretical point of view, study- 
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ing the quark-gluon plasma is a challenging task since it 
requires a deep understanding of QCD, and more gen- 
erally of gauge theories at finite temperatures. Several 
frameworks have been developed and have led to a great 
amount of works: Perturbative methods, potential mod- 
els, AdS/QCD duality, lattice QCD,. . . References about 
these topics can be found for example in the reviews . 

In principle, the most powerful technique to study non- 
perturbatively the properties of a quark-gluon plasma 
is lattice QCD. The equation of state of an SU(2) 
and SU(3) gluon plasma in lattice QCD were obtained 
in Refs. 0, 0, 0, [S, 9]. The equation of state of a 
quark-gluon plasma with non-vanishing fiavor number, 
Nf y^O, has also been computed in the more recent 
Refs. [13, UH [H, [13 ■ Notice that we focus here on the 
case where the chemical potential vanishes, but some re- 
sults have already been obtained at nonzero chemical po- 
tential (see for example Ref. [l3|)- In Fig.[Tl we show the 
equation of state obtained from pure glue SU(3) lattice 
computations in the continuum limit 7^ but in all cases 
{Nf — 0, 1, 2, 3), two important features are observed: (i) 
Energy and entropy increase sharply just after the phase 
transition temperature while the increase for the pressure 
is less pronounced, {ii) Energy and entropy seem to sat- 
urate below the Stefan-Boltzmann constant in the range 
T/Tc « 2 — 5. It has been argued in Ref. [lH| that finite- 
size effects were partly responsible of that behavior, but 
such numerical artifacts cannot explain the whole deficit. 
Actually, it is observed in the lattice computations of 
Ref. that the equation of state of the gluon plasma 
becomes compatible with the Stefan-Boltzman limit at 
very large temperatures. The ratio of pressure p/p^^ 
grows for example from 0.85 around Tc to a value compat- 
ible with 1 at T/Tc = 3 10^ [l^, that is an increase with 
a mean slope of order 10~^. Strictly speaking, the ther- 
modynamical quantities do thus not saturate below the 
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Stefan Boltzman constant. But, since we are interested 
in reproducing the lattice data of Fig. [1] for T/T^ < 5, 
and because the saturation rate is so small, we will fit 
the lattice data in the following as if that saturation was 
truly realized for T/Tc < 5. The error introduced by such 
an approximation is indeed completely negligible for our 
purpose. 




FIG. 1: (Color online) Energy density, entropy density, pres- 
sure, and interaction measure (or trace anomaly) of the gluon 
plasma versus T/Tc, as measured in pure glue lattice QCD at 
zero chemical potential Q (dashed lines) . The full horizontal 
line shows the Stefan-Boltzmann limit for a gas of massless 
transverse gluons. 



gluon mass) in order to reproduce the lattice data. Var- 
ious authors have proposed a procedure to treat statis- 
tically systems whose Hamiltonian depends on tempera- 
ture: In Ref. [2^, the authors start with the usual parti- 
tion function and simply replace the constant gluon mass, 
m, by a temperature-dependent one, m{T), which leads 
to an invariant expression for the pressure whereas the 
energy and entropy are modified in order to satisfy stan- 
dard thermodynamical relations between those quanti- 
ties. In Ref. [3l| however, the expression of the entropy 
is kept unchanged whereas both the energy and pres- 
sure are supplemented with an additional term involving 
dTm{T). Finally, in Ref. [Ill, the author proposes to 
keep the expression for the energy unchanged whereas the 
entropy and pressure get an additional term which also 
involves dTm(T). All these three formulations are obvi- 
ously not equivalent and there is still a debate to know 
which one, if any, is correct. In Sec. lIIH we propose a way 
to clarify the situation starting from the first principles 
of statistical mechanics while in Sec. llVl we show that all 
these formulations found in the literature demand sim- 
ilar temperature-dependent gluon masses to reproduce 
the lattice data. It implies that, at a qualitative level, 
those formulations are rather equally good. Moreover, 
we propose a new formulation where the expressions for 
the energy, entropy and pressure are invariant but where 
the Lagrange multiplier f3 is no longer equal to T^^ in 
order to ensure the so-called thermodynamic consistency, 
i.e. the fulfillment of the laws of thermodynamics. Note 
that we work in units where h = c = ks = ^, being 
Boltzmann's constant. 



B. Quasiparticle models 

There have been many attempts to understand the re- 
sults obtained in lattice QCD and to derive the equa- 
tion of state of a quark-gluon plasma from effective 
approaches. Indeed, QCD itself can be perturbatively 
solved only in the region of asymptotic freedom, i.e. for 
very high momenta or temperatures [T7| . But the con- 
vergence of the expansion in the strong coupling con- 
stant for the pressure is rather slow [l^ and conse- 
quently phenomenological models have been developed. 
There are mainly two frameworks: Strongly interact- 
ing quark-gluon plasma models taking explicitly into ac- 
count the pos sible existence of bound states beyond Tc 
EllSllil,!!!,!!!!!!,!!!], or quasipartlcle models, where 
the quark-gluon plasma is described as an ideal gas of 
massive bosons and fermioi is p , [26l . [27l [28l . ^29, 30, 31, 
m, m, m, HE H i3, il in this paper, we are 

mainly concerned with the quasiparticle formulation of 
a pure gluon plasma. As we recall in Sec. |TT1 it is nec- 
essary in such a framework to consider a phenomeno- 
logical temperature-dependent gluon mass (or thermal 



Most of the quasiparticle models reach the same con- 
clusion about the qualitative behavior of the thermal 
gluon mass: Just beyond the critical temperature Tc, 
the gluon mass has to be large and to decrease up to 
T/Tc ~ 1.5 — 2. Then, for even larger T, it increases 
essentially linearly. We interpret the large value of the 
gluon mass around Tc as a signal of strong color inter- 
actions among gluons. In Sec. |Vl we take those inter- 
actions into account by considering the gluon plasma as 
an ideal gas of glueballs and gluons, the color interac- 
tions being responsible for the formation of glueballs. We 
show that the lattice data can be reproduced if the ratio 
between the number of glueballs and gluons, n(T), de- 
creases monotonically when T increases so that a small 
value is reached when the temperature is larger than the 
dissociation temperature of the glueball. The glueball 
mass and dissociation temperature are computed using a 
spinless Salpeter equation (to take relativistic effects into 
account), the potential is obtained from lattice QCD [i^l 
and, to be consistent, the gluon mass is a linear function 
of T with the same slope as the one obtained from the 
asymptotic analysis performed in the quasiparticle for- 
malism. Finally, some conclusions and outlook are given 
in Sec. ED 
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II. IDEAL QUANTUM GAS OF BOSONS 
A. Constant mass 

As mentioned in the introduction, we study in the 
present work a pure gluon plasma in order to understand 
the available pure glue lattice data. In this section, we 
review briefly the reasons why an ideal gas of bosons with 
constant masses cannot reproduce lattice data, as well as 
why other attempts found in the literature, such as set- 
ting a momentum-cutoff for the gluons, are not suitable. 

From Fig. [H it is clear that the gluon plasma is not an 
ideal gas of massless bosons, that would only lead to the 
Stefan-Boltzmann constant. The first natural attempt to 
understand the equation of state obtained in lattice QCD 
is thus to consider that gluons have a constant nonzero 
mass. We give here the expressions of the energy density, 
Co, entropy density, sq, and pressure, po, in such a case 
since they will be useful later. For large enough volume, 
V, the sum over the possible quantum states is replaced 
by an integral and we have for a vanishing chemical po- 
tential (see for example Ref. (4ll . Chapter 5]) 



eo 



So 



Po 



dkk^q{e{k)lT) 



V 



X [kdke{k) + 3e{k)], (lb) 



Tso - eo 
d 



dk k\{e{k)/T)dke{k), 



(Ic) 



where d is the degeneracy factor, equal to 16 for trans- 
verse gluons (8 colors x 2 polarizations), and where q{x) 
is the Bose-Einstein distribution 



Using the following dispersion relation 

e(fc) = \/ k"^ + vn?, 
Eqs. ([1]) can be rewritten as 

eo(m, r)r"** 



(2) 
(3) 



—T dxq{x) \Jx'^ - A2, (4) 



so(m,r)r-3 ^ / dxq{x) v/a;2-v42 

67r2 J A 

X (4a;2 - A^) , (5) 

Po{m,T)T^^ = —T / dxq{x) (x^ - A 
67^2 J A 



(6) 



where 



The so-called interaction measure (or trace anomaly) is 
then found to be 

I{m,T)T-^ = (eo-3po)T-4 



d 
2^ 



dxq{x) \Jx'^ - A^. (8) 



Notice that / is positive (it is negative for non-relativistic 
bosons) and does not vanish as soon as the bosons have 
a nonzero mass. 
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Asymptote: 4.712 




3.344(r/r)"^ 




A{T) ^ m/T. 



(7) 



FIG. 2: (Color online) Energy density and interaction mea- 
sure of the gluon plasma versus T/Tc, as measured in pure 
glue lattice QCD at zero chemical potential [7]] (full circles 
and squares). The dashed-dotted lines outline the apparent 
saturation value at T/Tc < 5 and the behavior near Tc of 
the energy density. The solid line is a fit of the interaction 
measure beyond its maximal value (dashed line). 



This one-parameter - m - model cannot reproduce the 
lattice data even if there is some qualitative agreement. 
In particular, it cannot reproduce simultaneously some 
of the important characteristics summarized in Fig. [2] 

1. For T/Tc < 5, the energy density apparently sat- 
urates around 4.712, which is about 10% below 
the Stefan-Boltzmann constant dn'^/SO = 5.264 
(the gap is even of about 20% when one considers 
Nf ^ El), 

2. Around the critical temperature, the energy density 
increases very fast with a mean slope roughly equal 
to 26, 

3. The maximum of the interaction measure is located 
at T/Tc ~ 1.1 and its value is about 2.6, 

4. The decay of the interaction measure is given in 
good approximation by 3.344 {T/Tc)~'^ . 
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Indeed, a model with constant mass predicts that energy 
density, entropy density, and pressure will quickly sat- 
urate at the Stefan-Boltzmann constant. It is already 
enough to discard the model if we consider that the sat- 
uration occurs at the smaller value given above, namely 
4.712 instead of 5.264. However, even if we consider that 
lattice data did not yet saturate, the following arguments 
show that a model with a constant gluon mass is not able 
to reproduce these data since their four features enumer- 
ated above imply four very different values of the ratio 
m/Tc. At T/Tc = 4.5, the value of the energy dH) is 
equal to the lattice value 4.712 if A ~ 1, which implies 
m/Tc ~ 4.5. But the slope of the energy (|4]) can be 
bounded from above: 



< 



d 
2^ 
d 



^{^dAh{A)){-dT/TA) (9) 



m f T 

27r2^-^7; It; 



(10) 



Consequently, around Tc, the derivative of the energy 
density is smaller than around l.Sm/Tc, which implies 
m/Tc ~ 20 to reproduce the data. The interaction mea- 
sure ([8]) has a maximum located at A = 2.303 and the 
value of the maximum is equal to \.2d/{2ir^) ~ 1 (far 
from 2.6 obtained in lattice QCD). To reproduce the po- 
sition of the maximum obtained from lattice data we find 
that to/Tc ~ 2.5. At last, the decay of the interaction 
measure predicted by this model for large enough tem- 
perature is d/\2A^. Again to reproduce the data, we 
need m/Tc ~ 1.6. It is a simple exercise to verify that 
even with a degeneracy factor, d, considered as a free 
constant parameter, one cannot get a good quantitative 
agreement between this model and the lattice data in the 
whole available temperature range. 



B. Momentum cutoff 

Another version of this simple model is obtained by 
introducing a cutoff, K, for small momenta in Eqs. ([1]) 
as in Refs. [27l . [28j . The physical motivation of such 
a cutoff is that, near the critical temperature, gluons 
with low momenta should be bound into glueballs and 
should thus not contribute to thermodynamical quanti- 
ties related to an ideal gas of gluons. It is easy to see 
that this modified model predicts also a quick satura- 
tion of the energy density, entropy density, and pressure 
at the Stefan-Boltzmann constant. Moreover the addi- 
tional parameter K does not help to describe the large 
value of the derivative of the energy density around T^. 
This can be seen by using arguments similar to the ones 
given above when if = or by looking at the various fig- 
ures of Refs. [13, [ll]. Another possibihty is to consider 
a cutoff, if (T), which depends on the temperature as 
in Ref. Q. To have a relevant physical meaning, K{T) 
should be a decreasing function of T since one expects 
that, for high enough temperatures, glueballs will not 



be present in the plasma (see also Sec. |Vl. Using argu- 
ments similar to the ones developed to get the general 
form of m(r) in Sec. IIV Al one can prove that indeed, 
around Tc, wc have dTK{T) < 0. However, to reproduce 
the saturation of the thermodynamical quantities below 
the Stefan-Boltzmann constant, the cutoff must increase 
linearly with T for large enough temperatures. Qualita- 
tively, the shape of K{T) would be similar to the one of 
m{T), see Fig. [S] This behavior of K{T) for large T is 
problematic and its physical meaning is not obvious. 



C. Temperature-dependent mass 

For the reasons mentioned above, various authors have 
considered that a temperature-dependent gluon mass 
is the most relevant ingredient to be added to this 
model m, M, mi, [13, [11 [M M, [li- Indeed, since en- 
ergy density, entropy density, and pressure are decreasing 
functions of A, a simple way to make them saturate be- 
low the Stefan-Boltzmann constant for large T is to have 
a mass m[T) such that A — m{T)/T saturates to a non- 
vanishing constant in this regime of temperature. This 
implies that one must have 



m(T) - T for T > Tc 



(11) 



Although TO(r > Tc) is mostly a phenomenological pa- 
rameter that has to be fitted in order to reproduce the 
lattice results, the thermal gluon mass can be related to 
another important parameter characterizing a plasma, 
that is the plasma frequency. In a QED plasma for ex- 
ample, photons cannot propagate with a frequency be- 
low the plasma frequency. The situation is similar in the 
quark-gluon plasma: Gluons (plasmons) cannot propa- 
gate as free particles if their energy is too low. In fact, 
gluons acquire a thermal mass which is the plasma fre- 
quency. Perturbative calculations confirm that point: To 
leading order, the thermal gluon mass, also proportional 
to the Debye mass at large temperatures [33J, has been 
found to be proportional to y/ajT) T [H [il]. Fol- 
lowing standard notation, as{T) — g^(T)/47r with g{T) 
the strong coupling constant. Higher-order contributions 
have been calculated in the literature; it appears that the 
perturbative expansion converges very slowly and thus 
that the leading-order result is only valid at very high 



temperatures [4J] . Since as{T) ~ 1/ ^Jh\{T) varies 



more slowly than T, the linear behavior of the thermal 
gluon mass is found to be dominant at very large T in 
perturbative QCD. The constraint (fTT|) is thus in good 
qualitative agreement with already known results. 

It is worth mentioning that a resummation of the 
leading-order formulas leads to a modification of the ther- 
mal gluon mass that improves the convergence of the re- 
sults toward the Stefan-Boltzman limit. See in particu- 
lar Ref. 14511 ■ focusing on the convergence of the pressure, 
and Ref. [46j in which the hard thermal loop perturbation 
theory is successfully applied to reproduce lattice data at 
finite chemical potential within a quasiparticle approach. 
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Higher-order perturbative calculations also give good re- 
sults fH. 



III. STATISTICAL MECHANICS WITH 
TEMPERATURE-DEPENDENT HAMILTONIANS 

In a phenomenological model describing the gluon 
plasma, it is quite natural to assume that such a plasma 
can be seen as some gas of quasiparticles, or quasiglu- 
ons, as many authors have done in previous works (see 
the introduction). Standard statistical mechanics is the 
necessary framework to deal with gases of quasigluon. 
But, as we stressed in the previous section, a qualita- 
tive description of the pure glue lattice data demands 
the introduction of a temperature-dependent mass for 
the quasigluons. Such a thermal mass also emerges from 
QCD itself, as shown within perturbation theory. 

The problem is that standard statistical mechanics 
only deals with Hamiltonians which do not explicitly de- 
pend on temperature. There is a debate in the litera- 
ture about how it should be modified in such a case to 
be consistent. In this section, we propose an extension 
of statistical mechanics aiming at treating temperature- 
dependent Hamiltonians. Our procedure does not only 
allow to recover other existing formalisms in a unified 
way, it also leads to another new formulation which pre- 
serves the standard definition of energy, entropy, and 
pressure. We first consider classical systems in equilib- 
rium in the canonical ensemble and restrict our study 
to reversible processes to keep the discussion as simple 
as possible. Then the case of an ideal quantum gas of 
bosons, relevant for our study, is presented. 



A. General formalism 

Let us consider a probability density which is a func- 
tion of the Hamiltonian H{pi, qi,T): p = p{H). p is the 
density of probability of finding the dynamical variables 
of the system, namely {pi,qi\, within some volume of 
the phase space. We assume here that the Hamiltonian 
depends explicitly on the temperature T . The form of 
the function p can be determined using standard proce- 
dures; we give here the main ideas and refer the reader 
to standard textbooks for more details like, for example, 
Ref. [H, p. 50-55]. Consider two subsystems A and B 
with Hamiltonians Ha and Hb at equilibrium. The prob- 
abilities to find the subsystem A within the phase-space 
volume d\A and B within cIXb are given respectively by 
cIPa = Pa{Ha) d^A and by c?Pb = Pb{Hb) dXs following 
the definition of the probability density. Consider then 
the system C obtained by combining the two subsystems 
A and B. The probability to find both A within dXA and 
B within dXs , hence to find C within dXc = dXA dXs , 
is given by dPc = dPAdPs = pa{Ha) Pb[Hb) dXc- If 
the system C is also at equilibrium, we can write that 
dPc = PC {He) dXc- Neglecting the interaction between 



A and B, we have that He = Ha + Hb and 
Pc{Ha + Hb) = pa{Ha) Pb{Hb). 



(12) 



The intuitive idea underlying this constraint is that the 
description of a given system as a whole or as composed of 
two subsystems at equilibrium must be equivalent. Such 
a functional equation has a unique non trivial solution: 



Pa{x) = Cie pBix) 
and pcix) = CiCse-'^", 



-I3x 



(13) 



where Ci , C2 , and P are arbitrary constants with respect 
to the dynamical variables. Consequently, the general 
form of the normalized probability density is 



PiH) 



Je-f^HdX' 



m) 



(14) 
(15) 



where we have introduced the partition function Z al- 
though it does no longer play a central role in the present 
formulation. For the moment, /? is still arbitrary and will 
be fixed later. We recall that / dX denotes an integration 
on the phase space of the system. 

Another elegant way to obtain the expression ([M]) for 
the probability density is to use the so-called maximum- 
entropy estimate. According to Jaynes [i^, this is the 
least biased possible estimate on the available informa- 
tion. He showed that the form ((M]) maximizes the en- 
tropy of the probability distribution, i.e. the Shannon's 
measure [s^l, given by 



S = -In p 



plnpdX, 



(16) 



where we have introduced the notation x for the phase- 
space average of the quantity x. Notice that, in this 
formulation, /3 is a Lagrangian multiplier and, as above, 
does not depend on the dynamical variables. One advan- 
tage of this derivation is that it shows straightforwardly 
the intimate link between entropy and probability den- 
sity. The next quantity we consider is the energy, defined 
as the averaged Hamiltonian 



E = H = HpdX 



(17) 



The only unknown in the above relations is (3. We 
assume that this parameter depends only on the temper- 
ature T in the following way: 



T = 



1 



(18) 



The function f{x) in assumed to be positive, monotonic 
and such that f{x) S [0, oo[ so that the function /3(T) 
can be unambiguously defined. In order to determine 
/(/3), we use the laws of thermodynamics, that give a 
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relation between energy, entropy and /(/3). The first law 
links the variation of internal energy, of heat, Q, and 
of work, W, as follows: dE = SQ + 6W. For reversible 
processes we have 6W = —pdV. The second law re- 
lates the variation of heat with the variation of entropy: 
dS = 6Q/T (the equality holds for reversible processes 
only). Assuming that both the internal energy and the 
entropy are functions of /3 - thus implicitly of T - and of 
V, and combining the first and the second laws of ther- 
modynamics, we obtain 

dpEdp + dvEdV = T{d(}Sd(3 + dvSdV)-pdV. (19) 

Equating the terms in dP and dV, we find a relation 
between energy and entropy 



d0S^f{(3)dpE, 



(20) 



where we used Eq. (jlSp and, in the thermodynamical 
limit, the following expression for the pressure 



p = Ts - e = {TS-E)/V. 



(21) 

Equation ((20)) is the relation that allows to determine 
/(/3) and to know the relation between (3 and T through 
Eq. (fT5)) . Indeed, substituting expression into the 
definition p6|) of the entropy and using Eq. pT|) . we ob- 
tain 

5 = In Z + /?£;. (22) 
The derivation with respect to (3 of the entropy leads to 



E + fidfiE, 



(23) 



A comparison between Eqs. (|20p and ([23]) gives an equa- 
tion for /(/3): 



9gg(K,g„r=l//(/3)) 
dpE{T = 1/fm 



(24) 



This is, in general, a nonlinear first order differential 
equation for /. One thus gets f{/3, c) and, thanks to 
Eq. (I18|) . the function f3{T,c). The integration constant 
c can be constrained (sometime even fixed) by impos- 
ing that /(/3) is positive, bounded and monotonic for 
(3 e [0, oo[. Moreover, if there exists a temperature Tq 



such that dxH 



T=To 



= 0, the boundary condition needed 



to determine uniquely the solution of this equation is ob- 
tained by imposing that (3{To,c) = I/Tq in order to re- 
cover the usual formalism at this particular temperature. 
To clarify the procedure, we give an explicit example in 
Appendix [X] 

In this approach, the entropy, energy, and pressure are 
given respectively by their usual expressions ([TB)) . ([T7]), 
and (HJ). But, the dependence on T of the Hamilto- 
nian enforces a particular link between P and T, that 



can be found through the resolution of the nontrivial re- 
lation ([M)) . In standard problems, OtH = and one 
recovers the well-known link f{/3) ^ (3 ~ 1/T as a so- 
lution of Eq. (|24p . This general procedure has the seri- 
ous advantage of preserving the formal expressions of all 
the relevant thermodynamical quantities of the problem: 
The only modification arises at the level of the defini- 
tion of /3, which is not a physical parameter in itself. 
For computational applications however, this formalism 
is rather complicated since it needs an a priori knowledge 
of the solution of Eq. ([M]) if one wants to extract physi- 
cal informations about the system under study. Equation 
([24)1 can only be solved once the dependence on temper- 
ature of the Hamiltonian is explicitly known. However 
in the context of gluon plasma, the thermal mass of the 
gluons is unknown and must be determined from lattice 
data. That is why it is of interest to find more tractable 
ways of dealing with temperature-dependent Hamiltoni- 
ans. As we show in the following, our procedure allows to 
find such formulations, that corresponds to frameworks 
already in use in the literature. The study of the for- 
malism developed in this section will be the subject of a 
forthcoming paper. 



B. Alternative solutions 



First, notice that Eq. ([25)1 can be rewritten as 







dpS = dp{S+ udpH\p=,dv)^ l3dpE, (25) 

where we have introduced a new form, 5, for the entropy 
and where is some integration constant. The relation 
between the new entropy and the energy is then formally 
identical to the one given by Eq. (|20p provided we choose 
f{f3) = /3 as in standard statistical mechanics. The mod- 
ified pressure is given by Eq. (PTjl where S is replaced 
by S. In this formulation, the standard expressions for 
the energy and for (3 are preserved but the expressions for 
the entropy and pressure are modified. We then loose the 
usual connection between the probability density and the 
entropy. This formalism has been proposed in Ref. [s^]. 

Second, it is readily observed that another equivalent 
rewriting of Eq. ([23]) is 



dpS^PdpiE 



/3* 



di3H\i3=,dv) = l3di3E, (26) 



where we have introduced a new form, for the energy. 
Again, the relation between the new entropy and the en- 
ergy is formally identical to the one given by Eq. (PD|) 
provided that we choose also f{(3) = (3. The modified 
pressure is given by Eq. (PT|) where E is replaced by E. 
In this formulation, the standard expressions for the en- 
tropy and for (3 are preserved but the expressions for the 
energy and pressure are modified. In particular, the en- 
ergy is no longer the average of the Hamiltonian. This 



7 



formalism was first proposed in Ref . 13 1| and used in sev- 
eral other works, for example in Refs. [32l. [ssl. [3^. l46j . 

It is worth mentioning a third procedure that has been 
used in Ref. [2§| and consists in preserving the expres- 
sion of the pressure. We mention it for completeness but 
will not further study it in the present work. In can be 
deduced from Eq. ([2T|l that 



the energy density, entropy density, and pressure, i.e. 



T 
V 



InZ. 



(27) 



Consequently, one can leave the pressure invariant by 
setting /(/3) — (3 and by computing Z as usually done. 
But in this case, the entropy and the energy will be 
modified since the laws of thermodynamics demand that 
e = s/ (3 — p, with s = —P'^dpp: a term in dpH appears 
because of the dpp term. 

These three alternative solutions are derived from the 
laws of thermodynamics where the usual link f{/3) — (3 
is kept, but each one requires the standard form of the 
thermodynamical quantities to be modified. We think 
that the formalism developed in the previous section is 
the most fundamental one since it only demands a re- 
definition of /3, which is only a Lagrangian multiplier. 
Moreover, this redefinition of /3 as a function of T is only 
local since this is done through the differential equation 
m- ^/T = fm = (3 when dTH{p,,q,,T) = 0. In 
contrast, the corrections to the thermodynamic quanti- 
ties obtained from the alternative formulations derived in 
this section are non-local; they involve an integral over 
some range of temperatures, see Eqs. (f25l) and (|26|) . This 
implies that even if dTH{pi, g^, T) = in some large in- 
terval of temperatures, the corrections can be quite sig- 
nificant if the integration is performed over an interval 
of T for which the Hamiltonian depends on T. How- 
ever, the formalism proposed in Sec. IIII Al is far more 
complicated to deal with in numerical, phenomenolog- 
ical, applications when the dependence on temperature 
of the Hamiltonian is not known a priori. That is why the 
other approaches are useful too: From a computational 
point of view, it is easier to compute the extra integrals 
appearing in Eqs. (^5)1 or ([^ than to solve Eq. ([M)) if 
the temperature dependence of the Hamiltonian on T is 
not known. 



C. Ideal quantum gas of bosons 



The derivations presented in Sees IIII Al and [ITl Bl also 
holds formally in the quantum case provided that the av- 
erage takes correctly into account the statistics of bosons 
and fermions. We now focus on the case we are even- 
tually interested in: a gluon plasma. As we did in the 
classical case, we start with the standard expressions for 



eo 



So 



Po 



/ &fc2g(/3e(fc,/3))e(fc,/3), 

14 r dkeq{Pe{k,f3)) 
67r^ Jo 



(28a) 



X [A:9fce(fc,/3) + 3e(fc,/3)],(28b) 
/ dkk\{Pe{k,f3))dke{k,P), (28c) 



where e{k,i3) = e(fc, l//(/3)). We actually assume that 
the temperature dependence of e is known, the relation 
between T and P being given by Eq. (|18|) which was ob- 
tained from the laws of thermodynamics. 

It is convenient to write the entropy as follows 



So = InZo + (3eo, 



where 



lnZo = -— r/ dkk^\n{l-e 
27r^ Jo 



,/3£(fc,/3) 



The derivative of sq with respect to /? leads to 

where the average is now given by 
d 



dkk^q{l3e{k,(3))dpe{k,l3). 



(29) 



(30) 



(31) 



(32) 



The comparison between Eq. (PT|) and Eq. (|20p leads to 
the following expression for /(/3): 



dpe{T = 1/fm 
dpeo{T = l//(/3)) 



(33) 



This formula for / is formally identical to the one ob- 
tained in the classical case; only the definition of the aver- 
age is different. The same comments about the boundary 
condition can be made. 

It is also possible to recover in the quantum case 
the other formalisms proposed in the literature. In the 
formulation where the expression of the energy density 
is preserved, the expression for the new entropy den- 
sity is still given by Eq. ([^5]) but the average is now 
defined by Eq. ((32)) and the pressure can be obtained 
through Eq. (HI]). A similar remark applies for the for- 
malism which preserves the expression of the entropy, see 
Eq. (|26p . Consequently, what we call Model 1 in the next 
sections is defined by the equations 



e(i)=eo, s(i) =so + S(i\ p = po + 



13 - 



with 



(34a) 



(34b) 
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Model 2 is similarly defined by 



= (2) _ 



60-^(2), s = so, p^po + B'-'), (35a) 



(2) 



with 



l3 



B^'>iP) = / dpe\f3=,diy. 



(35b) 



In both Models 1 and 2 we have (3 = l/T as usual, and 
Co, sq: and po are given by Eqs. ((28)) . The integration 



constants, [3'^^ and are arbitrary parameters. The 
function e{k,l3) is given by 



;(fc, /?) = e(fc, l/T) = ^k^+m^{T), 



(36) 



with m{T) the thermal gluon mass. If that function was 
unambiguously known, Models 1 and 2 would lead to in- 
equivalent results and the general solution ([55)1 should 
rather be used. However, m(T) is a parameter of the 
model, that can be fitted on some lattice data. Its be- 
havior at large T should nevertheless be coherent with 
Eq. PT|) . in agreement with perturbative QCD. Thus it 
can be expected that all the presented formalisms should 
lead to very similar results provided that the gluon mass 
and the integration constants, /3i , are properly chosen. 
In particular. Models 1 and 2 are compared in Sec. IIVI 
where we show that they lead to similar temperature- 
dependent gluon masses and that they are both able to 
reproduce the lattice data with excellent accuracy. Con- 
sequently, from a practical point of view, they are both 
rather equivalent because there is not much constraints 
on this thermal gluon mass. 



IV. APPLICATION TO THE GLUON PLASMA 
A. General features of m{T) 

In this section, we give some general characteristics of 
the function m{T). Let us consider the Model 1 where the 
expression of the energy is conserved (similar conclusion 
can be obtained with Model 2 by using the expression 
of the entropy). From the lattice data, we know that 
the slope of the energy near the critical temperature is 
positive and rather large (see Fig. [2]) . From Eq. (|4]) we 
find 

Bt/tS^^T-^) - ^{dAh{A)){dT/TA) > 0. (37) 

Since dAh{A) < 0, we find that ra{T)/T is a decreasing 
function of the temperature 



m(T) 



< 0. 



(38) 



From Eq. (U, it is easy to see that in order to reproduce 
the saturation of the energy for T/T^ > 4, the thermal 
mass m(T) should be a linear function of T: m(T) = 



in T, i.e. A = rh, with m = 0.973. We can also show 
that the large mean slope of the energy near the critical 
temperature implies that the derivative of m(T) should 
be negative in this region. Indeed, 

c't/t, (eoT""*) < ma.x{-dAh{A))ma,x{-dT/T,A). 

(39) 

From the expression ^ of h{A), it is readily computed 
that max(— 9a/i(j4)) = M ~ 1.6. Let us assume that 
dTm{T) > everywhere, then 



iiiax.{—dj'/j'^A) — max 



< max 



< max 



1? 








m{T) 








m{T) 









T 



5tto(T) 



(40) 



where rh{T) = mT for T > T and rh{T) = mT for T < T 
and where T is the temperature where the thermal mass 



should be linear (T/Tc - 2 



< 



' 3). Consequently we obtain 



-Mm- 

2^2 T,, 



T 



< 1-3— < 4, 



(41) 



since T/Tc ~ 1, T/Tc ~ 2 - 3, m ~ 1 and M ~ 1.6. This 
upper bound on the derivative of the energy is much lower 
than the value given by the lattice data. Large values for 
this derivative can only be obtained if dTm{T) < near 
the critical temperature. 

In consequence, we have analytically shown that an 
agreement with lattice QCD can be obtained provided 
that 

dTm(T > Tc) < 0, and m{T > = 0.973 T. (42) 



B. Numerical results 

In the previous section, we have shown that the shape 
of m{T) is rather constrained by the lattice data within 
the frameworks of Models 1 and 2. Let us now explicitly 
extract numerically m(T) from these data. 

We begin by considering Model 1, where the form of 
the energy density is preserved and consequently given 
by Eq. (|28ap in which the dispersion relation is cho- 
sen. At a given temperature, T* , the thermal gluon mass, 
m{T*), can be obtained by numerically solving the equa- 
tion e(i'(T*,m(T*)) = eo(T*,m(r*)) = e^'^^T*), where 
e'^*(T*) is the lattice energy density at the considered 
temperature. The computed gluon mass is plotted in 
Fig. [3] and can be well fitted by the following form 



m(i)(r) 

T, 



T 



mi 



[T/T, - 



(43a) 
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with 



form 



(43b) 



mo = 0.873, mi = 0.612, 
m2 = 0.983, ma = 0.411. 



The observation of Fig. [3] and of the fitted form ((43)) 
clearly shows the different behaviors predicted in 
Sec. lIV Al First, the linear increase of m(T) is obvious for 
T/Tc > 2.5, and corresponds to the region III in Fig. [31 
However, the slope mo differs from the asymptotic value 
of 0.973 predicted in the previous section by about 10%. 
This can be understood by remarking that ms is rather 
small while mi is of the same order of magnitude than 
mo: The term supplementing the linear one in Eq. (|43ap 
still brings a non negligible contribution at large temper- 
ature, causing the fitted slope mo to be smaller than in 
the case of a genuine linearly rising mass, see Eq. 
Second, m{T) strongly decreases for T/Tc ~ 1.0 — 1.2, 
corresponding to region I in Fig. [31 The fitted form 
we get is actually singular near the critical temperature, 
the parameter playing the role of a critical exponent. 
Third, there exists an intermediate zone between the sin- 
gular and the linear behaviors, in which m(T) reaches a 
minimum. This zone corresponds to region II in Fig. [31 

In Model 2, the form of the energy density is no longer 
preserved as in Model 1, but the entropy density is so. 
That is why the thermal gluon mass at a given tem- 
perature, T* , can be numerically computed in a very 
similar way by solving the equation s'^^^T* , m(T*)) = 
soiT*,m{T*)) = s'^*(r*). We recall that sq is given by 
Eq. (|28bp . The computed thermal gluon mass is plotted 
in Fig. [31 and can be accurately fitted by the following 



5.0 
4.5 
4.0 



P 3.5 



Model 1 
Model 2 
Model 1 - Fit 
Model 2 - Fit 




FIG. 3: (Color online) Thermal gluon masses obtained by 
fitting Model 1 (circles) and Model 2 (triangles) to the lat- 
tice data of Ref. [3| , see Fig. [1] Models 1 and 2 are defined 
by Eqs. (|34|) and (|35p respectively, with the dispersion rela- 
tion p6)) . The fitted forms (|43)) and (|44l) are also plotted 
(solid lines). 



with 



m(2)(r) T ki 



ka = 0.724, fci = 0.982, 
fca = 0.973, fcg = 0.345. 



(44a) 



(44b) 



Equation ([44]) is formally equivalent to ([43]); only the 
values of the numerical coefficients are slightly different. 
Thus the same comments as for Model 1 can be done. It 
is worth noting that the regions where the linear increase 
and strong decrease occur are identical within Models 1 
and 2. A physical interpretation of the thermal gluon 
mass can thus be given independently of the considered 
Model. 
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FIG. 4: (Color online) Interaction measure of the gluon 
plasma computed with Model 1 (solid line) and Model 2 
(dashed line), and compared to the lattice data of Ref. 0] 
(circles). Model 1 is defined by Eqs. (|34p with the gluon 
mass ((43} and Model 2 is defined by Eqs. (|35|l with the gluon 
mass (|44p . The dispersion relation ((36} and the values ((45} 
for the integration constants are used. 

The thermal gluon mass has been fitted on one of 
the three thermodynamical quantities available in lattice 
QCD: energy density for Model 1, and entropy density for 
Model 2. The remaining quantities can now be numer- 
ically computed within both Models by using Eqs. ((3^ 
and ((35|) with the dispersion relation ([35]). provided that 
the integration constants ensuring an optimal agreement 
with lattice QCD are known. The following fitted values 



0.435, and Tc/3?^= 0.445 



(45) 



lead to an excellent agreement with the available lattice 
data, as it can be seen in Figs. (Hand O 

The interaction measures computed with Model 1 and 
Model 2 are quasi indistinguishable from each other and 
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FIG. 5: (Color online) Same as Fig.[T] but the results obtained 
with Model 1 are also plotted for comparison. Model 1 is 
defined by Eqs. (|34|l with the gluon mass (|43|l . the dispersion 
relation (|36|l . and the value (|45p for the integration constant. 



Color interactions above Tc 



We have considered up to now that the gluon plasma is 
an ideal boson gas, where gluons are transverse and free 
but have a temperature-dependent mass m(T). From 
Sec. IIV Al we can conclude that reproducing the lattice 
QCD results demands first that m(T > Tc) T and 
second that m{T > Tc) decreases fast enough (see Fig.[3|). 
Such a nontrivial behavior can be intuitively explained by 
invoking color interactions above T^. It is indeed widely 
accepted that, at the critical temperature Tc, the medium 
undergoes a phase transition and becomes deconfined. It 
does not mean however that the color interactions vanish: 
They are actually screened because of the great amount 
of color charges in the medium, and the residual potential 
is no longer confining as for T = 0. These residual color 
interactions can be rather important for T/T^ — 1 ~ 2, 
as suggested by several lattice QCD studies [4^ . 

One can think about the mean field approximation to 
have a first guess about the influence of screened color 
interactions. In this picture, the gluon dispersion relation 
should be modified as follows: 



from the lattice data. Again, both formalisms lead to 
nearly identical results. That is why, for clarity, we have 
only plotted the results of Model 1 in Fig. [5l The curves 
computed with Model 2 would have been indistinguish- 
able from those of Model 1. Notice that in our approach, 
the thermal gluon mass is fitted so that the asymptotic 
behavior of the interaction measure corresponds to lat- 
tice QCD, that is e - 3p cx (see Fig. [2]). Such a 
quadratic increase is compatible with previous theoret- 
ical results 15 ll and with the more recent unquenched 
lattice study |52l |. It is worth mentioning that other ap- 
proaches rather favor e — Bp cx T [ssl . [Hj, Issl ] . The in- 
teraction measure is thus a quantity that would deserve 
further studies since there is not yet a general agreement 
concerning its asymptotic growing. 

The values we find for the integration constants are 
almost equal: Their average value is Tc(3 — 0.44, corre- 
sponding to the temperature T/Tc ~ 2.27. This is the 
typical temperature at which rrS^\T) = m'-^-' (T), as it 
can be seen in Fig. [31 It is not a coincidence: Mod- 
els 1 and 2 are designed to reproduce the same data. 
Then if to(i)(T) = to(2)(T), the only way for both Mod- 
els to give identical results is to have B^^\p) = 0, thus 
/Si^-* — /Sp-* = P, in rough agreement with the fitted val- 
ues (PS)) . The integration constant is thus not really a free 
parameter since its value can constrained by the thermal 
gluon mass once Model 1 and Model 2 are compared. 

Finally, it is important to stress that the terms B'^'^\ 
given by Eqs. p4b|) and ((35b|) . are not small corrections 
as one could have thought. Without these terms, even 
the qualitative behavior of the various thermodynamical 
quantities is wrong. 



v//c2 + m2(T) = Vfc2 + m2(T) + V{T), (46) 

where V{T) is the effective mean potential energy felt 
by a gluon, and where fn{T) is a priori different of 
m{T). When T becomes very large, it can reasonably 
be assumed that V{T) vanishes. Consequently, one has 
to(T ^ Tc) = m{T 3> Tc) — rfiT. Our main physi- 
cal assumption is then the following: Since fri(T) is the 
thermal gluon mass in a temperature range when the 
gluons are free, it can be seen as the rest mass of a 
free gluon for any T > T^. Let us now consider that 
y{T > Tc) > m(T > Tc), i.e. that the color interactions 
become dominant near Tc. Then, by squaring Eq. (|46p . 
one gets m^{T) = m^{T) + V^{T) + 2y'k^ + fh^{T)V{T) 
and consequently m2(T > Tc) « V'^{T > Tc) since the 
potential term dominates the right hand side. 

We are thus led to the following interpretation for 
m(T): At large T (region III in Fig. [3]), to(T) tends to 
the rest mass of a free gluon in the gluon plasma accord- 
ing to the perturbative QCD result. But near Tc (region 
I in Fig. [3]), the behavior of m{T) is dominated by the 
existence of non negligible screened color interactions. 
Region II in Fig. [3] is finally a transition regime in which 
these interactions progressively vanish. According to this 
scenario, screened color interactions play an important 
role for T/Tc ~ 1 — 2, and one can wonder whether glue- 
balls can form or not at these temperatures. The next 
section is devoted to answer to this question. Notice that 
the possible glueball formation above Tc has already been 
suggested in Ref. [svt as a mechanism explaining the sud- 
den increase of the effective gluonic degrees of freedom 
near Tc that is observed in this last work. 
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V. EXISTENCE OF GLUEBALLS ABOVE Tc 
A. Effective Hamiltonian for glueballs 

In a constituent gluon (or quasiparticle) picture such as 
the one we develop here, a glueball is a bound state of at 
least two gluons. Let us focus on two-gluon glueballs. Be- 
ing the lightest and presumably the most strongly bound 
ones, they should be the easiest glueballs to be produced 
in the gluon plasma. In a deconfined medium, a binary 
gluon state may exist in several colored configurations fol- 
lowing the decomposition 8(8)8 = 1®8®8©10©T0®27. 
As the strength of color interactions is proportional to the 
color Casimir operator of the gluon pair, the last three 
configurations are irrelevant as far as glueball formation 
is concerned since they lead to interactions which are 
either vanishing (10, 10) or repulsive 27 [iO, HI]- How- 
ever, both the singlet and octet configurations lead to 
attractive interactions, these interactions in the singlet 
channel being twice as large as in the octet one. The 
most favored glueball from an energetic point of view is 
thus a two-gluon bound state with the gluon pair in a 
color singlet. The dynamics of the gluon pair also comes 
into play at this stage: The most strongly bound gluon 
pairs will be those with a minimal value of the radial 
quantum number (n = 0) and of the orbital angular mo- 
mentum. If the gluons were longitudinal the minimal 
value of the square orbital angular momentum would be 
{L"^) = for the O"*""*" state. However, we have seen that 
the large-T behavior of the gluon plasma is compatible 
with transverse gluons. In this case, (i^) = 2 is the min- 
imal allowed value, corresponding to the 0*+ glueballs as 
shown in Ref . [57,] . 

Denoting the static potential between a color-singlet 
quark- antiquark pair by y(r, T), a relativistic Hamilto- 
nian describing the aforementioned lightest glueballs is 
the following spinless Salpeter one 

2v/p^ ■ — ^ 



G 



rn?{T) 



'V{r,T), 



where = + (^L-^) jr^ and where 



m T) T T 
——^ = Th—^ 0.973 — 

Tr T, Tc 



(47) 



(48) 



is the free gluon mass introduced in the previous sec- 
tion. The 9/4 factor comes from the color Casimir op- 
erator. Such a Casimir scaling for the static energy be- 
tween sources in various color re pre sentation has been 
confirmed by the lattice study [5^. We choose for 
m the value that reproduces that saturation value of 
the thermodynamical quantities following the analysis of 
Sec. lTOI Notice that, since Tc is estimated to be around 
270 MeV by pure glue lattice calculations Q], one has 



m(T) = 0.263— GeV. 



(49) 



We point out that building an effective glueball Hamil- 
tonian by starting from a best-known quark-antiquark 



one has already led to a successful description of glue- 
balls at T = f57|. That is why we find relevant to 
apply it in this case also. A last remark has to be done: 
Our framework leads by construction to the same mass 
for the scalar and pseudoscalar glueballs. At T = this 
degeneracy can be lifted by the introduction of instanton- 
induced forces \5T\. Such forces are not taken into ac- 
count by the present model and one can expect that the 
ground state of Hamiltonian (|47|) . whose mass is denoted 
Mc{T), is rather an average mass of the scalar and pseu- 
doscalar glueballs. Although the current understanding 
of this topic is far from being complete, we can neverthe- 
less mention that instantons effects might be less impor- 
tant at high temperatures following Ref. [s^. Actually, 
we have checked that the results that we obtain in the fol- 
lowing do not demand an accurate knowledge of Mq{T). 

A key ingredient in Hamiltonian (|47)) is the potential 
energy V{r, T) between a static quark-antiquark pair. It 
is well known from lattice QCD that this potential is 
compatible with a funnel shape ar ~ b/r dX T — [60j . 
but the situation is less clear when T > 0. The potential 
energy that is the most readily obtained in lattice QCD 
is the quark-antiquark free energy F{r, T) [1, [6l[. We re- 
call that, thermodynamically speaking, the free energy of 
a system is the energy that is available in the system to 
produce a work once the energy losses due to the increase 
of the entropy have been subtracted. As also noticed in 
Ref. (20| . in a potential approach however, the potential 
energy of the system should be the total energy that it 
contains, no matter it will be lost or not in heat trans- 
fers. Such a potential energy corresponds to the internal 
energy of the system, usually denoted hy U = F + TS, 
where S is the entropy. The internal energy of a quark- 
antiquark pair is thus the quantity we choose as potential 
term. It has been computed in lattice QCD in Refs. [40| : 
we give a plot of these results in Fig. [6] Notice that 
those Nf = computations are the most relevant for our 
purpose since we consider a genuine gluon plasma. 

It can be checked in Figs. [5] and [7] that the lattice data 
are accurately fitted by the following form 



Vir,T)^-aiT) e"'''^) + c(r). 



where 
a{T) 

c{T) 



ao 



Co 



02 In 



and 



T 



(50a) 



T 



,6(T) = 60 + 61-, 



(50b) 



ao = 0.459 GeV, ai = 0.915, aa = 1.159 GeV, 

60 = 0.111 GeV, 61 = 0.489 GeV, 

Co = 0.341 GeV, ci = 0.808. (50c) 

It has to be stressed that the form ((5(7| is the one that 
gives the best fit of the lattice data but is not moti- 
vated by any physical theory predicting such a form. We 



12 




T/Tc = 1.13 
T/Tc = 1.18 
T/Tc = 1.25 
T/Tc = 1.40 
T/Tc = 1.65 
T/Tc = 1.95 
r/Tc= 2.60 
T/Tc = 4.50 
T/Tc = 7.50 
- Fits 



4 6 
r (GeV"') 



10 



FIG. 6: (Color online) Internal energy of a static quark- 
antiquark pair computed in lattice QCD for different values 
of T/Tc and for Nf = (dots). Lattice data are taken from 
Refs. S'lid compared to the fitted form (|5Up for some values 
of T/Tc (solid lines). 
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FIG. 7: (Color online) Values of a(T), b{T), and c(T) obtained 
by a fit of the lattice QCD data to the form (|50al) (dots), 
compared to the analytical curves (|50b[l (solid lines). 



just use it as a convenient parameterization of the lat- 
tice QCD results. It clearly appears from Eq. (|50ap that 
the potential energy is no longer confining. Exponential 
potentials indeed only admit a finite number of bound 
states. 



B. Numerical results 

All the terms appearing in Hamiltonian (j47p are now 
explicitly known, and its ground state mass can be nu- 
merically computed. To this aim we use the Lagrange 
mesh method, which is a numerical procedure allowing 
in particular to accurately solve eigenequations associ- 
ated to relativistic Hamiltonians [6^. The evolution of 
the lowest-lying glueball mass with the temperature is 
given in Fig. [S] The numerically computed evolution of 
the glueball mass with T is accurately fitted by the form 



Mg{T) = -c{T)+2m{T)+2h{T) e{T) for T < 1.13 T^, 



where 



£(r) 



£o = 0.818, 



sqT/T, - ei 



0.921, 



0.958. 



(51a) 
(51b) 
(51c) 



The glueball mass we find is around 1.8 GeV at T « Tc, 
then increases to reach a maximal value of about 2.8 GeV. 
Notice that the mass near the critical temperature is sim- 
ilar to the one obtained at zero temperature [1^ . To our 
knowledge, the behavior of the scalar glueball mass ver- 
sus the temperature has not been studied a lot in the 
literature. We can nevertheless quote the lattice study of 
Ref. [64| that finds a reduction of 20% of the scalar glue- 
ball mass when one goes from T — to T = Tc, and the 
more recent work |65| finding an almost constant glue- 
ball mass from T = to T = Tc. Beyond the qualitative 
behavior of Mq{T), an important result we find is that 
the ground state is bound up to T = 1.13 Tc and then 
dissociates in the medium above this temperature. Nu- 
merically, the dissociation temperature is reached when 
the binding energy of the system vanishes. Our model 
thus predicts the existence of glueballs in the tempera- 
ture range T/Tc = 1 ^ 1.13, but the existence of bound 
states is a very stringent criterion: Glueball resonances 
can indeed appear in the continuum even if the gluons 
are not bound. Following the lattice results of Ref. [1^1 , 
glueball resonances can even be expected up to 1.9 Tc- 

Since glueballs can be present in the deconfined 
medium, we propose to recompute the thermodynami- 
cal properties of the gluon plasma by assuming that it is 
a mixing between an ideal gas of transverse gluons and 
an ideal gas of glueballs; the glueball abundance n{T) 
depending on the temperature. This last approach will 
be referred to as Model 3; it shares with Model 1 the 
property that (3 — 1/T and that the form of the en- 
ergy is preserved. Its spirit is a bit similar to the one 
of the hadronic resonance gas model, assuming that the 
hot hadronic medium can be described as an ideal gas 



made of all possible resonance species (see Refs. [66l |67| 
for more informations). We actually consider that the 
screened color interactions "generate" color singlet scalar 
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FIG. 8: (Color online) Numerically computed lowest-lying 
glueball mass, that is the ground state mass of Hamilto- 
nian (gZl, versus T/Tc (solid line). The fitted form ((51]) is 
also plotted for comparison (dashed line). The curve stops at 
the glueball dissociation temperature, namely 1.13 Tc- 

and pseudoscalar glueballs in a first stage, and that these 
glueballs behave as free particles in the gluon plasma in 
a second stage. Consequently, if 

eo[d,m,f3) = — ^^-^ — dk (52) 

^'^ Jfim e — i 

is the energy density of an ideal gas of bosons with mass 
771 and with d degrees of freedom, then the total energy 
density of the mixed gluon-glueball gas is 

e(3) = [l-77(T)]eo(16,?7iT,l/r) 

+7i(T)eo(2,AfG(T),l/r), (53a) 

where two degrees of freedom are associated to the glue- 
ball gas, accounting for the lowest-lying 0*"'" states. A 
priori, n(T) should vanish above 1.13 Tc because glue- 
balls are then not bound anymore. However, two-gluon 
resonances can in principle appear in the continuum 
above the dissociation temperature. The simplest way 
to take this phenomenon into account is to allow n{T) 
to be nonzero above the dissociation temperature. In 
this sector, formula (|51ap remains well-defined, roughly 
simulating a gluon pair in the continuum. 

The unknown function 7i(T) can be computed by fit- 
ting Eq. (|53ap to the lattice energy. The result is given in 
Fig. [HI it appears that the numerically computed curve 
is accurately described by the following form 

„(r) = e-«o(T/T.-i)"i^ (53b) 

with no — 3.358 and ni = 0.541. The glueball abun- 
dance is nearly 100% a,t T — Tc, then decreases to reach 
33% at 1.13 Tc, the dissociation temperature of the two- 
gluon glueballs. But as we said previously, resonances are 
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FIG. 9: (Color online) Glueball abundance computed by fit- 
ting Eq. ((53a)) to lattice QCD versus T/Tc (solid hue). The 
fitted form (|53b[) is plotted for comparison (dashed line). 

then still expected to form in the continuum, justifying a 
nonzero glueball abundance at higher temperatures. Fi- 
nally, 77(T) is less than 5% at 1.9 Tc- Such a negligible 
value is coherent with the fact that glueball resonances 
are expected to disappear above that temperature [65| . 
We have checked that the quantitative behavior of n{T) 
is not very sensitive to the glueball mass, Ma{T). The 
key result of Model 3 is rather that cq (16, 7tiT, 1 /T) alone 
is unable to fit the available data and consequently that 
an additional term accounting for glueballs is needed. 

The entropy density can be computed from Eqs. ((53aP 
and (gni)- It reads 

/•l/Tc 

s(3) ^ / I3dpe^^\l/I3)dl3. (53c) 

Jl/T 

The upper bound of this last integral ensures that 
s(3)(Tc) = 0, in qualitative agreement with lattice QCD. 
Finally, the pressure can be computed thanks to the def- 
inition ((2T|). that is 

p(3) ^r5(3)_g(3)^ (53d) 

The results are plotted in Figs. (TU) and (TTl and compared 
to lattice QCD. As it was the case for Models 1 and 2, 
Model 3 leads to an excellent agreement with the lattice 
data, although relying on a different physical picture of 
the gluon plasma. 

VI. CONCLUSIONS AND OUTLOOK 

It is now widely accepted that the equation of state of 
the gluon plasma, coming from pure gauge lattice QCD 
computations, can be accurately reproduced by model- 
ing the gluon plasma as a gas of transverse gluons with 
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FIG. 10: Same as Fig. [l] but lattice data are this time com- 
pared to Model 3 (solid gray line) defined by Eqs. (|53p . 




FIG. 11: (Color online) Same as Fig. [5] but lattice data are 
this time compared to Model 3 (solid gray line) defined by 
Eqs. ([53ll. 



a temperature-dependent mass. As we have outlined in 
the beginning of this paper indeed, such a quasiparti- 
cle model is indeed in disagreement with lattice QCD 
if a constant gluon mass is used. One is thus led to 
deal with temperature-dependent Hamiltonians. In that 
case, standard formulas in statistical mechanics have to 
be modified in order to enforce the thermodynamical con- 
sistency, but the procedure to achieve such a task varies 
from one work to another. In the frameworks that can 
be found in the literature so far, the standard expression 
of only one thermodynamical quantity can be preserved 
in order to enforce the thermodynamic consistency, i.e. 
to satisfy the laws of thermodynamics. The expressions 
of the other quantities have to be modified: Either the 
pressure 0|, the entropy [3l|, or the energy [3^ is kept 



invariant. 

In this work, we have clarified the situation by show- 
ing that all the existing formulations can be derived in 
a simple unified way. In the process, we have uncovered 
a new possible formulation for which the standard form 
of each thermodynamical quantity is preserved but for 
which (3 is no longer equal to 1/T. The function P{T) 
has to be extracted from a first order nonlinear differ- 
ential equation expressing the fulfillment of the laws of 
thermodynamics. We think that this last formalism is the 
most fundamental one, since it only demands a change in 
the definition of the Lagrangian multiplier /3, which has 
no physical meaning a priori. Moreover, the corrections 
to standard statistical mechanics implied by this new for- 
mulation are only local in T - i. e. they vanish in regions 
where the Hamiltonian does not depend on T - while cor- 
rections found in other formulations are non-local in T. 
However, this new formulation is far more complicated to 
deal with in numerical applications when the dependence 
of the Hamiltonian on temperature is not known. That 
is why the other approaches are also useful to study the 
quark-gluon plasma. 

Consequently, we focused on two formulations: The 
ones that preserve the form of the energy and of the 
entropy. It can be analytically shown that, indepen- 
dently of the considered formulation, reproducing the lat- 
tice data leads to constraints on the thermal gluon mass, 
m{T). It must be strongly decreasing just after the crit- 
ical temperature and grow linearly asymptotically. A 
numerical fit of the thermal mass on the available data 
confirms this behavior and eventually leads to an excel- 
lent agreement with lattice QCD. Both frameworks lead 
to nearly indistinguishable results as expected, and to 
very similar thermal gluon mass. 

Mean-field-inspired arguments show that the singular 
behavior of the thermal gluon mass near Tc accounts 
for residual color interactions, which are still strong in 
the early stages after deconfinement. The potential en- 
ergy coming from such screened color interactions has 
already been computed in lattice QCD, allowing us to 
build a consistent Hamiltonian describing the interac- 
tions between two transverse gluons in a color singlet, 
that is the channel in which the color interactions are the 
strongest. It appears that the two-gluon ground state, 
corresponding to the scalar and pseudoscalar glueballs, 
remains bound up to T = 1.13 Tc- We have then pro- 
posed a last description of the gluon plasma, in which 
this medium is seen as a ideal mixture of free gluons and 
colorless glueballs. The agreement with lattice QCD is 
as good as with the previous approaches, with a glueball 
abundance that is very large near the critical tempera- 
ture, takes the lower value of 33% at the dissociation tem- 
perature of the hghtest glueballs, and becomes negligible 
after 1.9 Tc, where even continuum glueball resonances 
are expected to disappear (65j . This interpretation of 
the gluon plasma draws a bridge between the quasiparti- 
cle approach and other models focusing on the existence 
of bound states after deconfinement [20| . 
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From an experimental point of view, the main result of 
the present study is the prediction that the gluon plasma, 
and thus presumably the quark-gluon plasma, might be a 
glueball-rich medium in the early stages after deconfine- 
ment. This brings support to previous studies arguing 
that an important amount of glueballs can be formed 
in relativistic heavy ion collisions [H, [6^ . The experi- 
mental detection of the scalar glueball in the quark-gluon 
plasma could be achieved through the scenario developed 
in Refs. [T^I which roughly suggests that, although the 
bare scalar glueball would be nearly stable in the quark- 
gluon plasma, it should mix with scalar mesons. Then, 
such a "physical" glueball, denoted as G, could decay 
mostly in the channels G — > tttt and G ^ 77 through its 
mesonic component, leading to an enhancement of the 
number of events versus the two-photon (or two-pion) 
invariant mass. In our model, the bare glueball mass is 
mostly located around 2.8 GeV; a peak in the 77 or tttt 
channels can thus be expected not too far of 2.8 GeV, de- 
pending on the strength of the meson-glueball coupling. 

We finally stress that, if quarks were included in our 
model, the number of bound states above the critical tem- 
perature would increase since mesons, diquarks, quark- 
gluon states, etc. can also form. We leave the extension 
of our approach to the full quark-gluon plasma for future 
works. The effects of a nonzero chemical potential will 
also be leaved for subsequent studies. 
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APPENDIX A: DETERMINATION OF /3(T) FOR 
TEMPERATURE-DEPENDENT 
HAMILTONIANS: AN EXAMPLE 

In order to illustrate the general procedure given in 
Sec, nil Al we study the particular case of a classical ideal 
gas with temperature-dependent mass m{T), or equiva- 
lently 'm{f{j3)) because of the definition (fT8| . The Hamil- 
tonian reads 



H = 



2m(/(/3))' 

and one finds that, for system of N particles, 
^ r27rm(/(/3))T^^/' 



V 



(Al) 



(A2) 



where V is the volume of the system. The normalized 
probability density is thus known and it can be computed 



that 
E 



3N 
2^' 



and dpH — 



' 2/3 m(/(/3)) ■ 



(A3) 



where the prime denotes a partial derivation with respect 
to the argument of the considered function. These last 
two equalities allow to rewrite Eq. as 



^ rn {f{/3)) , 



0. 



(A4) 



Let us consider the following form for m{T) to illus- 
trate the procedure 



(l^^l+4A7T)/2 



m{T) = rriQe 



with (5 > 0. The equation for /(/3) then reads 



(A5) 



SP^fiP) = (/(/?) - /3)Vl + W(/?). (A6) 

To determine uniquely the function /, we need a bound- 
ary condition. For T — )■ 00, the mass tends to mp which 
is constant. In this limit, one recover the standard sta- 
tistical mechanics and /(/3) = (3 = 1/T. Consequently 
the boundary condition is /(O) = 0. The unique solution 
of the nonlinear differential equation (|A6p is then 



(A7) 



The relation between f3 and T is then given by 



-1 + ^1 + 45/T 



25 



(AS) 



For T 45 (or (5 ^ 0), we just recover the standard 
relation (3 — 1/T. In this formalism, we also find that 



E 
N 



35 



-1 + ^1 + 45/T 



(A9) 



We can now compare this last energy formula with the 
energy obtained within Models 1 and 2. For the Model 
1, where the expression for the energy is preserved, we 
simply have the standard expression E/N = 3T/2 while 
for the Model 2, where the expression for the entropy is 
preserved, the energy takes the form (remember that in 
this formalism /3 = 1/T) 



E 
N 



'ST 


/■ 


~2 ^ 


J/3* 


3T 


3_f 


T + 




3T 


36 


T + 


2 



dpH\i3=^du, 

^ d,m{v) 1 ^ 
—aiy, 



(AlO) 



In 



1 + VI + ^5v 



-1 + Vl + 



The correction to the energy {E/N — 3T/2) can be com- 
pared for each formalism. Of course, for Model 1, this 
correction is vanishing; in this case corrections would be 
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FIG. 12: (Color online) Comparison between the corrections 
to the energy for the new formalism and Model 2 as a function 
of T/5. We used T*/5 = 1. The evolution of m{T)/mo is also 
presented. 



associated to the entropy. Consequently, in Fig. [121 we 
compare only corrections to the energy obtained with 
Model 2 and with the new formalism proposed in this 
paper together with the evolution of the mass as func- 
tions of the temperature T/6. 

We notice that the corrections to the energy from 
Model 2 (and corrections to the entropy from Model 1) 
are non-local since they involve integrals over some range 
in temperature, see Eqs. (^5)) and (^5)1 . This means that 
those corrections are still significative in regions where 
the mass is essentially constant (in this example the cor- 
rections are logarithmic in T) while the corrections to the 
energy from the new formalism are essentially localized 
around the region where the mass depends significantly 
on the temperature. This is indeed what we expect: If the 
Hamiltonian does not essentially depend on T over some 
large interval of temperature, the statistical mechanics in 
this interval of T should be essentially the same than the 
standard statistical mechanics. 



[1] J. C. Collins and M. J. Perry, Phys. Rev. Lett. 34, 1353 
(1975). 

[2] E. V. Shuryak, Phys. Rep. 61, 71 (1980). 

[3] BRAHMS Experiment, Nucl. Phys. A 757, 1 (2005) 
nucl-ex/0410020 ; PHOBOS Experiment, Nucl. Phys. 
A 757, 28 (2005) nucl-ex/0410022 ; STAR Experi- 
ment, Nucl. Phys. A 757, 102 (2005) nucl-ex/0501009 ; 
PHENIX Experiment, Nucl. Phys. A 757, 184 (2005) 
[nucl-ex/0 410003 . 

[4] U. Kraemmer and A. Rebhan, Rept. Prog. Phys. 67, 
351 (2004) hep-ph/0310337; H. Satz, J. Phys. G 32, 
R25 (2006) hep-ph/05122ir ; E. Shuryak, J. Phys. G 
35, 104044 (2008) ar Xiv:0804.1373j C. DeTar, PoS 
LAT2008, 001 (2008) [arXiv:0811.2429 . and references 
therein. 

[5] J. Engels, F. Karsch and H Satz, Nucl. Phys. B 205, 545 
(1982). 

[6] J. Engels, J. Fingberg, K. Redlich, H. Satz and M. Weber, 

Z. Phys. C 42, 341 (1989). 
[7] G. Boyd, J. Engels, F. Karsch, E. Laermann, C. Lege- 

land, M. Liitgemeier and B. Petersson, Phys. Rev. Lett. 

75, 4169 (1995) |hep-lat/9506025 . 
[8] G. Boyd et a l, Nucl. Phys. B 469, 419 (1996) 

hep-lat/9602007 . 
[9] M. Okamoto et ai, Phys. Rev. D 60, 094510 (1999) 

hep-lat/9905005 . 
[10] F. Karsch, E. Laermann and A. Peikert, Phys. Lett. B 

478, 447 (2000) hep-lat/0002003 . 
[11] A. AU Khan et al, Phys. Rev. D 64, 074510 (2001) 

hep-lat/0103028 . 
[12] Y. Aoki, Z. Fodor, S. D. Katz, and K. K. Szabo, JHEP 

0601, 089 (2006) h ep-lat/05 10084 . 
[13] C. Bernard et al, Phys. Rev. D 75, 094505 (2007) 

hep-lat/0 611031| ]. 
[14] C. Bernard et al., Phys. Rev. D 77, 014503 (2008) 

arXiv:0710.1330 . 
[15] F. Gliozzi, J. Phys. A 40, F375 (2007) [hep-Iat/070i020] . 



[16] G. Endrodi, Z. Fodor, S. D. Ka tz, and K. K. Szabo, PoS 

LAT2007, 228 (2007) [arXiv:0 710.4197 . 
[17] H. D. Politzer, Phys. Rep. 14, 129 (1974). 
[18] K. Kajantie, M. Laine, K. Rummukainen, and 

Y. Schroder, Phys Rev. D 67, 105008 (2003) 
hep-ph/0211321 ]. 
[19] V. M. Bannur, Phys. Lett. B 362, 7 (1995). 
[20] E. V. Shuryak and I. Zahed, Phys. Rev. D 70, 054507 

(2004) hep-ph/0403127 . 
[21] J. Liao and E. V. Shuryak, Nucl. Phys. A 775, 224 (2006) 

hep-ph/0508035 . 
[22] J. Liao and E. V. Shuryak, Phys. Rev. D 73, 014509 

(2006) hep-ph/0510110 . 
[23] B. A. Gelman, E. V. Shuryak and I. Zahed, Phys. Rev. 

C 74, 044908 (2006) nucl-th/0601029 . 
[24] B. A. Gelman, E. V. Shuryak and I. Zahed, Phys. Rev. 

C 74, 044909 (2006) nucl-th/0605046 . 
[25] J. Liao and E. Shuryak, Phys. Rev. C 75, 054907 (2007) 

hep-ph/0611131 . 
[26] D. H. Rischke, M.I. Gorenstein, H. Stocker and W. 

Greiner, Phys. Lett. B 237, 153 (1990). 
[27] D. H. Rischke, M. I. Gorenstein, A. Schafer, H. Stocker 

and W. Greiner, Phys. Lett. B 278, 19 (1992). 
[28] D. H. Rischke, J. Schaffner, M. I. Gorenstein, A. Schafer, 

H. Stocker and W. Greiner, Z. Phys. C 56, 325 (1992). 
[29] V. Goloviznin and H. Satz, Z. Phys. C 57, 671 (1993). 
[30] A. Peshier, B. Kampfer, O. P. Pavlenko and G. Sofi, 

Phys. Lett. B 337, 235 (1994). 
[31] M. I. Gorenstein and S. N. Yang, Phys. Rev. D 52, 5206 

(1995). 

[32] A. Peshier, B. Kampfer, O. P. Pavlenko and G. Sofi, 

Phys. Rev. D 54, 2399 (1996). 
[33] P. Levai and U. Heinz, Phys. Rev. C 57, 1879 (1998) 

hep-ph/9710463 . 
[34] R. A. Schneider and W. Weise, Phys. Rev. C 64, 055201 

(2001) hep-ph/0105242 . 
[35] D. Zwanziger, Phys. Rev. Lett. 94, 182301 (2005) 



17 



hep-ph/0407103 . 

P. Castorina and D. Zappala, Europhys. Lett. 80, 22001 
(2007) hep-ph/0703218; . 

P. Castorina and M. Mannarelli, Phys. Rev. C 75, 054901 
(2007) hep-ph/0701206 . 

V. M. Bannur, Phys. Lett. B 647, 271 (2007) 
liep-ph/0608232 ; Eur. Phys. J. C 50, 629 (2007) 
jhep-ph/0508069 ; Phys. Rev. C 75, 044905 
[2007) l iep-ph/0609188 ; JHEP 0709, 046 (2007) 

(hep^ph70604158 . ^ 

V." Chandra and V. Ravishankar, arXiv:0812.1430 
O. Kaczmarek, F. Karsch, P. Petreczky, and F. 
Zantow, NucL Phys. Proc. SuppL 129, 560 (2004) 
■hep-lat/ 0309121 ; P. Petreczky, Eur. Phys. J. C 43, 51 

(2005) hep-lat/0502008 . 

L. D. Landau and E. M. Lifshitz, Statistical Physics 

(Pergamon Press, Bristol, 1958). 

E. Shuryak, ZhETF 74, 408 (1978). 

J. Kapusta, Nucl. Phys. B 148, 461 (1979). 

K. Kajantie et al, Phys. Rev. Lett. 79, 3130 (1997) 

;hep-ph/9708207. 

J. P. Blaizot, E. lancu, and A. Rebhan, Phys. Rev. D 68, 
025011 (2003) hep-ph/0303045'. 

A. Rebhan and P. Romatschke, Phys. Rev. D 68, 025022 
(2003) hep-ph/0304294 . 

M. Laine and Y. Schroder, Phys. Rev. D 73, 085009 

(2006) hep-ph/0603048 . 

J. D. Walecka, Fundamentals of Statiscal Mechanics (Im- 
perial College Press and World Scientific Publishing Co. 
Pte. Ltd., 2000). 

E. T. Jaynes, Phys. Rev. 106, 620 (1957). 

C. E. Shannon, BeU System Tech. J. 27, 379 (1948); tbid., 
623 (1948). 

R. D. Pisarski, Prog. Theor. Phys. Suppl. 168, 276 (2007) 
hep-ph/06 12191| . 

M. Cheng et al, Phys. Rev. D 77, 014511 (2008) 
arXiv:0710.0354 . 

D. Zwanziger, Phys. Rev. Lett. 94, 182301 (2005) 



Phys. Rept. 443, 55 (2007) 



hep-ph/0407r03 
[54] D. E. Miller, 

hep-ph/0608234 . 
[55] F. Giacosa and R. Hofmann, Phys. Rev. D 76, 085022 

(2007) hep-th/0703127 . 
[56] N. Boulanger, F. Buisseret, V. Mathieu and C. Semay, 

Eur. Phys. J. A 38, 317 (2008) arXiv:0806.3174 . 
[57] V. Mathieu, F. Buisseret and C. Semay, Phys. Rev. D 

77, 114022 (2008) ar Xiv:0802 .0088 . 
[58] S. Gupta, K. Huebner and O. Kaczmarek, Phys. Rev. D 

77, 034503 (2008) arXiv:0711.2251 . 
[59] T. Schafer and E. V. Shuryak, Rev.^ Mod. Phys. 70, 323 

(1998) hep-ph/9610451 . 
[60] G. S. Bah, Phys. Rept. 343, 1 (2001) [hep-ph/00013T2] . 
[61] L. D. McLerran and B. Svetitsky, Phys. Rev. D 24, 450 

(1981). 

[62] C. Semay, D. Baye, M. Hesse, and B. Silvestre-Brac, 

Phys. Rev. E 64, 016703 (2001). 
[63] C. J. Morningstar and M. J. Peardon, Phys. Rev. D 60, 

034509 (1999) hep-lat/9901004 ; Y. Chen et al, Phys. 

Rev. D 73, 014516 (2006) hep-lat/0510074 . 
[64] N. Ishii, H. Suganuma and H. Matsufuru, Phys. Rev 

66, 094506 (2002) hep-lat/0206020j. 
[65] X. F. Meng et al, arXiv:0903.1991 
[66] E. Beth and G. E. Uhlenbeck, Physica 4, 915 (1937); 

Hagedorn, Nuovo Cim. A 56, 1027 (1968). 
[67] F. Karsch, K. Redlich and A. Tawfik, Eur. Phys. J. 

29, 549 (2003) hep-ph/03 03108, ; Phys. Lett. B 571, 

(2003) hep-ph/0306208 . 
[68] C. Gao, X. Li and W. Lu, Phys. Rev. C 52, 421 (1995). 
[69] I. N. Mishustin, L. M. Satarov and W. Greiner, J. Phys. 

G 32, L59 (2006) hep-ph/0606251 . 
[70] y. Vento, P hys. Rev. D 75, 055012 (2007) 

[hep-ph/060"92r9 ; V. Mathieu, N. Kochelev and 

V. Ve nto, Int. J. Mod. Phys. E 18, 1 (2009) 

larXiv:0810.4453 . 



D 



R. 

C 
67 



